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Abstract 
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and the universality class is different from the Z]\f invariant Villain model. At large 
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1 Introduction 



Recently [15| we have analysed in detail the deconfining phase transition in the SU{2) 
Georgi-Glashow (GG) model in 2+1 dimensions. The mechanism of confinement in this 
model at zero temperature is due to the "plasma" of the monopole-instantons and is 
well understood The model is weakly interacting all the way up to the deconfining 
temperature, which allowed us to study the phase transition quantitatively. We found 
that taking into account the excitations of the heavy charged particles was crucial for the 
correct description of the transition. The transition is associated with the restoration of the 
magnetic Z2 symmetry ^ in accordance with general arguments of |Q . The universality 
class of the transition was found to be 2d Ising. 

Whereas for SU{2) gauge theory there is overwhelming consensus that the transition 
should be in the universality class of the Ising model, the situation is much less clear for 
large A^. The point is that for > 3 one can write down different 2d spin models, and they 
have different critical behaviour. For example the A^-state Potts models have first order 
phase transition for > 4 while Villain models have second order transition which is 
of the BKT type, and is thus in the universality class of U{1) Whether the transition 
in the SU{N) gauge theory is similar to either one of those, is an open interesting question. 

In this paper we consider a general Georgi-Glashow type SU (N) gauge theory, where 
at zero temperature the gauge group is spontaneously broken to f/^~^(l). Just like the 
SU{2) GG model, the theory is weakly interacting. At zero temperature it is confining, 
and the monopole "plasma" description of confinement has long been known 0. It has 
also been studied from the point of view of magnetic Zn symmetry in [§. 

Our main finding is that the transition in the model is second order, and is distinct 
from that of Villain model. Although we are unable to identify the fixed point theory with 
a known two dimensional conformal theory, we argue that the relevant model at large A^ 
must be a deformation of a theory with a large value of the UV central charge c = 0{N), 
which may be SU{N)i WZNW model. 

The paper is structured as follows. In Section 2 we describe the model as well as 
the monopole and magnetic symmetry based approaches to its low energy dynamics. In 



Section 3 we derive the dimensionally reduced model relevant for the study of the phase 
transition, and discuss the role of the heavy charged particles. In Section 4 we study the 
transition with the help of the renormalization group analysis in the reduced theory. We 
show that the RG equations have a self dual infrared fixed point. We explain why the 
GG model close to the transition does not behave like Villain model, even in the range of 
parameters where one might expect it to do so. In Section 5 we point out to similarities 
between the behaviour of some quantities in the GG model close to criticality and in the 
hot Yang Mills theory. Finally in Section 6 we discuss our results. 

2 The model 

We consider the SU{N) gauge theory with scalar fields in the adjoint representation in 
2+1 dimensions. 

C = -^tiF^^F'"' + trL>^$L>^$ - V{^) (1) 

where 

= A^T- F^, = d^A, - d,A^ + g[A^, A,] 
$ = $«T'^ D^^ = d^^ + g[A^,^. (2) 

are traceless hermitian generators of the SU{N) algebra normalised as tr(T"T'') = 
Depending on the form of the Higgs potential, there can be different patterns of 
gauge symmetry breaking. Since most of the details of the potential are unimportant for 
our purposes, we will not specify it except for restricting it to the region of the parameter 
space where classically the gauge symmetry is broken to the maximal torus 

SU{N) ^ C/(l)^-i (3) 

We also restrict ourselves to weakly coupled regime, which means that the ratios M-^jg^ 
are large for all N'^ — N massive VF-bosons. 
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2.1 The perturbative spectrum. 

To characterise the perturbative spectrum of the theory it is convenient to use the Cartan- 
Weyl basis {H^,E"'), where iJ* generate the Cartan subalgebra which is of the dimension 
of rank of ^f/( AT): r = N - 1. 

[H\H^] = z,j e[l,2,..N-l] (4) 

and E"' are the N{N — 1) ladder operators which satisfy 

[H\E'^]= a'E'^, (5) 
[E", E^] = N^j E'^+^ if a + /3 is a root (6) 
= 2a-H ii a = -(3 (7) 

The A^ — 1 dimensional root vectors a = (a^, a^, ...a^~^) form the dual Cartan subalgebra. 
There are obviously N{N — 1) such vectors corresponding to dim{SU{N)) —Ta.nk{SU{N)) 
but only — 1 of them are linearly independent. The non- vanishing inner products in the 
Cartan- Weyl basis read as 

tr(if\ H^) = ^6'\ tr(E", E^) = ^S'^'-^. (8) 

At the classical level A^ — 1 gauge group generators are unbroken, which we choose to 
correspond to {H^)- Therefore classically there are A^ — 1 massless photons and A^(A^ — 1) 
charged massive W-bosons. 

Our Weyl basis is chosen in such a way that the Higgs VEV is diagonal. Since the 
matrix $ is traceless, there are A^ — 1 independent eigenvalues. In terms of the A^ — 1 
dimensional vector h = (/ii, /;,2, /is, ../iAr-i)0 we have 

<^>=h-H, A^ = A^-H + Y.^tE^ (9) 

d 

^ For concreteness we order these numbers hi > h2 > .../iAr_2 > hp^^i, which also breaks the discrete 
Weyl group. 
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For concreteness let us choose the following basis for the Cartan subalgebra; 

H, = ^diag(l, -1, 0, ...0), H2 = ^diag(l, 1, -2, 0...0) 

... giv-i= , ^ diag(l,l,l,...l,-(iV-l)) (10) 
V2iV(iV - 1) 

As long as /i ■ a 7^ for all roots, the gauge symmetry is maximally broken. The masses of 
the W-bosons can be read off from the second term in the lagrangian 

ghrlA,, $ ]^ = ^ E Af^Afh^h.aW (11) 

=^Ms = g\h-d\ (12) 

The W-bosons corresponding to the A^ — 1 simple roots i = 1, N — 1 (arbitrarily 
chosen set of linearly independent roots) can be thought of as fundamental, in the sense 
that the quantum numbers and the masses of all other W-bosons are obtained as linear 
combinations of those of the fundamental W-bosons. These charges and masses are 

Q0 = g0, M^ = gh-^ (13) 

As an example consider the case of S'?7(3) broken down to U{1) x f/(l). There are 6 massive 
W-bosons. The simple roots can be taken as 

The remaining non-simple positive root is 

«3 = /3i-/32 = (1,0). (15) 

The other three roots are — —a^. The masses of corresponding W-bosons are 

Mw, = |(/^i + V3h2), Mw, = |(/ii - V3h2), Mw, = gh (16) 

for hi > v^^2- Observe that if /i2 = 0, two of the masses become degenerate. In this 
case SU{3) is still broken down to U{1) x f/(l) since all three masses are non-vanishing 
but the spectrum is invariant under an additional Z2 symmetry. This Z2 symmetry is the 
charge conjugation with respect to the charge H2, which interchanges the roots (3i and /?2- 



In general though, this charge conjugation symmetry is broken by the VEV of Higgs |[TT[| . 
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2.2 The monopole-instantons and the Polyakov effective La- 
grangian. 

Non-perturbatively the most important contributions in the theory are due to the 
monopole-instantons. Those are classical, stable, finite action solutions of the Euclidean 
equations of motion arising due to the nontrivial nature of the second homotopy group of 
the vacuum manifold {Il2{SU (N) /U{1)^~^) = Z^). The magnetic field of such a monopole 
is long range. 

— * 

B^ = ^.3-H (17) 
The — 1 dimensional vectors g are determined by the non-Abelian generalisation of the 



Dirac quantisation condition p3|, 12 



e'^^-" = I (18) 
Solution of this quantisation condition are 

g = -Y.n.f3: (19) 

where /3* are the dual roots defined by /3* = We will be working with roots 

normalised to unity, and thus (3* = (3. The integers rij are elements of the group 112 O- 
The monopoles which have the smallest action correspond to roots taken once. The action 
of these monopoles in the BPS limit is 

M,= l^A.a-=l^ (20) 

9 9 

Just like with ly-bosons we can think of monopoles corresponding to simple roots as 

fundamental ones with magnetic charges and action 

47r ^ 47r ^ ^ 

Qi = —Pi Mi = —h- A (21) 
9 9 



For example, in the case of SU{3) (see eqs. dl^JTq) ) the monopole action spectrum (in the 
BPS limit) is 

Stt Stt 47r 

Ml = —{hi + V3h2), M2 = —{hi - V3h2), M3 = —hi. (22) 

9 9 9 
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The effect of these monopoles is to impart finite mass to all the perturbatively massless 
"photons" . The derivation of the effective Lagrangian follows exactly the same lines as the 
original derivation of Polyakov for the SU{2) theory The resulting low energy effective 
theory is written in terms of the — 1 component field, rf, with the following Lagrangian 

^eff = ^,{d,vf + E ^exp(^a . ff) (23) 

The sum is over all A^(A^ — 1) non vanishing roots. The potential induced by the monopoles 
is proportional to the monopole fugacity 

Ml = ^„ = constant 24 

9^ 9 



e(^) is such that 1 < e < 1.787 and e(oo) = 1. 



The photons at weak coupling are obviously much lighter than the W - bosons and 
thus are the only relevant degrees of freedom in the low energy sector. 

2.3 The magnetic symmetry. 

The global symmetry structure is very important for the understanding of the deconfining 
transition. The relevant symmetry in the present model is the magnetic Zn symmetry. 
We now wish to explain how this symmetry is implemented in the effective low energy 
Lagrangian. Our discussion parallels the SU (2) case [|15| . 



The order parameter of the magnetic symmetry is the set of magnetic vortex operators 
Vi, i = 1, ...,N — 1. These operators were constructed explicitly in P]. These operators 
carry magnetic fluxes of the A^ — 1 U{1) Abelian magnetic fields. The defining commutation 
relation for Vi is 

All 

[V,{x), B{y)] = tS,V{x)6\x - y) (25) 

9 

Here B is the — 1 dimensional vector of magnetic fieldsQ, whose j-th component is the 

projection of the non-Abelian field strength onto the direction of the Cartan subalgebra 

^Note that these magnetic fields can be constructed in gauge invariant way from the non-Abehan field 
strengths and the Higgs field, see |3|. 



generator Hj, and Wj are — 1 weight vectors of SU (N). The choice of the — 1 out of 
weight vectors is arbitrary. Change in this choice will lead to the redefinition of the vortex 
operators such that the new operators will be products of the old ones and their conjugates. 
It is always possible to choose these weights so that together with the "fundamental" roots 
Pi they satisfy the relation 

Wipj = ^Sij (26) 

The flux eigenvalues in eq. ( ]25|) are dictated by the requirement of the locality of the vortex 
operators and is analogous to the Dirac quantisation condition. The explicit form of the 
vortex operators in terms of the field i] in eq.(p3D is 



V,{x) = ^e^^^ (27) 

with 

Xi = w^-f]^f] = 2j2 f^iXi (28) 

i 

The effective Lagrangian can be written as a nonlinear cr-model in terms of Vi as 



with A ^ oo. The matrix A^j = 2f3* ■ (3j depends on the choice of the fundamental roots. 
With the conventional choice of positive roots, where (3il3j = —1/2, i ^ j, it is the Cartan 
matrix of the Lie algebra. All its diagonal elements are equal to 2, while all its off diagonal 
elements equal to —1. We will find it however more convenient in the following to use a 
different set of fundamental roots, for which = 1/2, i ^ j- Such a choice is always 
possible for any SU{N). For this choice of roots the off-diagonal matrix elements of Aij 
are all equal to 1. 

For SU{3) we have 

^=(12) PO) 
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and the effective Lagrangian 

HiiYiV; + c.c) + UV1V2 + c.c) + ^3(^1 ^2' + c.c) (31) 

The magnetic symmetry has an obvious and simple representation in this effective 
Lagrangian as — > exp{27rm/A^}l^. 

As long as only small ffuctuations of the phase fields Xi are important, the Lagrangian 
eq.(|3l[) is equivalent to the eq.(|^). Thus at low temperature the descriptions based on 
these Lagrangians are equivalent. The difference appears only when the phase nature of 
Xi plays a role. Indeed, since Xi are treated in eq.(pTD as phases, dynamically one allows 
configurations in which these phases have nontrivial winding. On the other hand in eq. ( p3D 
such configurations cost infinite amount of energy. As discussed in detail in and |TB| 
the winding configurations correspond to the heavy ly-bosons. In fact the explicit relation 
between the vorticity of the fields Vi and the electric charges is given by 

-g"^'^ ^^'^'-^'^ ^^^^ 



+00 



Thus the difference between the two Lagrangians is important whenever the physics of the 
W bosons plays an important role. We have seen in the case of the SU (2) theory that W^s 
are indeed important near the phase transition temperature. The same turns out to be 
true for arbitrary SU{N). We thus have to be careful to treat the W - bosons properly in 
the transition region. In the next section we will set up this treatment. 

3 The reduced theory. 

Throughout this paper we are working in the weak coupling regime and thus the photon 
masses in eq.(^) are exponentially small. Thus already at very low temperature (T oc M^) 
one can use the dimensionally reduced version of the theory, since all the thermal modes are 
significantly heavier than the zero Matsubara frequency mode. Since the critical temper- 
ature for the deconfining transition is of order g"^ (see [|^), we can safely use dimensional 
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reduction close to the transition. The zero Matsubara frequency sector is described by the 
two dimensional Lagrangian 

^-ff = ^M^^^^^)' + T^^^P(*« ■ ^) (33) 

However, as we noted before, our description should include W bosons, and so the fields 
T] should be treated here as phases with periodicity appropriate to eq.(^). In fact the 
Lagrangian also has to be augmented by a four derivative "Skyrme" term, which fixes 
the energy of the winding states to be equal to the masses of W bosons ||15|. We can 



however simplify things further, by noting that the density of W bosons at criticality is 
exponentially small due to the Boltzmann factor suppression. Thus W^s can be treated in 



the dilute gas approximation in the same way as was done in To do this explicitly we 
first have to understand how to write partition function in the presence of one W boson of 
a particular type. 

Let us first consider a W boson corresponding to one of the fundamental roots 
Using eq.(|l3D, eg. (^6]) and eg. (|32|) we see that this W boson corresponds to unit vorticity 



of the field Vk and zero vorticity of all other fields Vj, j ^ k. To create such a vortex in 
the path integral we must introduce an external "current" which forces the discontinuity of 
the field Xk 

Xk = Xk + 27r (34) 
The partition function in the presence of one W boson is thus 
Z = I D[xix)] 

exp |- y" d^y ^Q^2rp - J^iy, x)){d^Xj - Jf,{y, x)) + Y,Ca cos(z2 " ■ AXi) 

with 

Jliy, x) = 27i5ikn^{y)5iy E C^) (35) 

with Cx a curve that starts at the location of the vortex (the point x), and goes to infinity, 
and is the unit normal to this curve. The insertion of this current forces the normal 
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derivative of Xk to diverge on curve C, so that Xk jumps by 2tt across C. Since in the 
rest of the space Xk is smooth, the path integral is dominated by a configuration with unit 
vorticity of Xfc0- 

The path integral eq.(|35[) differs from the partition function in the vacuum sector by 
the linear term in the Lagrangian 

- J '^^y i^i ■ i^i d^XiJf, = dx^e^^ d^0k ■ V (36) 

Defining in the standard way the dual field fj, 

idpiV = ^t^udyT] (37) 

we can recast the contribution of this particular W boson in the form of the following extra 
term in the Lagrangian 

- ^^^^ ■ ^ (38) 

This procedure can be repeated for W boson corresponding to an arbitrary root a with 
the only difference that in eg. (P^) , the root [3k is replaced by the root a. To create several 
ly-bosons one just inserts the external current which is the sum of the currents creating 
individual VF's. 

Dilute ensemble of such objects with small fugacities /Xq, is then given by 

^ = n E / n n dy^^^^ (39) 



The summation over the number of W^s can be easily performed, see |T^. The result is 
the partition function with the Lagrangian 

^eff = g^^2y (^M^)^ + E Caexp(m " rf) + ^ /i«exp(i^5 ■ fj) (40) 

with summation in both terms going over all non- vanishing roots of SU{N). The coeffi- 
cients proportional to the fugacities of the corresponding W bosons 

fia(xexp{-MwjT} (41) 



^Note that even though expHcitly depends on the curve Cx, the partition function itself does not, 
since changing the integration variable Xii^) ~* Xii^) + 27r, x d S where the boundary of S* is — 
is equivalent to changing C — x into in the definition of the current. 
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Eq.(^D]) is the dimensionally reduced theory which we will now use to study the phase 
transition. 



4 The phase transition. 

4.1 Monopoles versus charges. 

To study the phase transition we may first attempt to disregard the W boson induced term 
in the effective Lagrangian. If we do that, we are back to the theory eq.(^). This theory is 
easily analysed. The first interesting thing about it is that since the group is simply laced 
(all the roots are of unit length) the anomalous dimensions of all the interaction terms are 
equal. The scaling dimension of all the monopole induced terms is 

AttT 

Am = (42) 
This immediately tells us that at the temperature 

Tbkt = ^ (43) 

ZTT 

all these interactions become irrelevant. Thus at Tbkt one expects the Berezinsky- 
Kosterlitz-Thouless transition to take place. Above this temperature the infrared behaviour 
of the theory is that of — 1 free massless particles. Note that Tbkt does not depend on 
the number of colours A^. If the picture just described where true, the universality class of 
the phase transition would be that of f/^^^(l). 

This of course is exactly the same situation as encountered in [|l4l in the SU (2) case. 
Again just like in SU (2) case this conclusion is incorrect due to the contribution of the W 
bosons. To see this it is simplest to ask what would happen at high temperature if there 
were no monopole contributions at all. This amounts to studying eq.(0) with = 0. 
This theory describes non-compact electrodynamics with A^ — 1 photons and the spectrum 
of charged particles given by eq.(0). This limit is again simple to understand, since the 
theory is exactly dual to the theory with monopoles and without charges. The scaling 
dimensions of all the W induced perturbations are equal and are given by 

A. = £^ (44) 
12 



Thus the perturbations are irrelevant at low temperature, but become relevant at 

T„c = |^ (45) 

Since Tjyc < Tbkt this tells us that we can not neglect the effects of charges at criticality. 
The story of SU (2) exactly repeats itself. Even the value of the temperature at which the 
scaling dimensions of the charge- and monopole induced perturbations are equal does not 
depend on A^. 

We expect therefore that the actual transition temperature is 

Tc = 1^ (46) 

at which point all perturbations have the same scaling dimension. This expectation is 
confirmed by the renormalization group analysis. 

4.2 Renormalization group analysis 

The renormalization group equations for the theory eq.(50) were studied in [16]. In general 



the equations are quite complicated due to the cross correlations between different opera- 
tors. For this reason the space of parameters of the theory has to be enlarged if one wants 
to study the flow whose UV initial condition is provided by eq.(^) with arbitrary values 
of fugacities. However there is one simple case, that is when the initial condition is such 
that all the monopole fugacities are equal = = ^, and all the charge fugacities are 
equal /io,. = /Iq,^ = fi. This initial condition is stable under the RG flow. On this subspace 
the RG equations, written in terms of the scaled temperature t = and dimensionless 
fugacities, read 

^ = 2n'Nt{„' - e) (47) 
1^ = (2 - 1)/. - 27r(iV - 2)^^' (48) 
^ = {2-t)C- 2niN - 2)e (49) 
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These equations have exactly the property reflecting our previous discussion. That is the 
points t = 2, /i = and t = 1/2, ^ = are both unstable. The stable IR fixed point is 

to = l f^o = Co= ^^^^_^^ (50) 

One can in fact easily check that in the three dimensional space of couplings t, ^ and this 
point has two attractive and one repulsive direction. This is precisely what one expects 
from the IR fixed point located on the critical surface, the two attractive directions being 
the tangential directions to the surface. 

The RG equations have an obvious duality symmetry, /i — ^, t ^ 1/t. This is the 
refiection of the transformation rj fj on the level of the Lagrangian eq.(^). The points 
t = 1, /i = ^ are symmetric under duality, and this ensures existence of a self dual fixed 
point. This is important, since the exact position of the fixed point is scheme dependent. 
Its existence however is assured by the duality symmetry. 



What is the nature of this fixed point? For = 2 we were able in |T^ to fermionize 
the fixed point theory and show explicitly that it is equivalent to one massless Majorana 
fermion. We are not able to perform a similar analysis for arbitrary A^. There are however 
several comments that we would like to make. Phase transitions in Zj^ invariant spin 
models have been studied quite extensively. A nice recent discussion of the situation is 
given in |T7|. One considers a spin model of one phase field 6 with a symmetry breaking 



term of the type h cos{N6} which breaks the U{1) symmetry down to Z^. When the 
coefficient h of this symmetry breaking term is large, the model resembles Potts model and 
thus (for A^ > 4) has a first order phase transition. When the breaking is small on the other 
hand, the behaviour is similar to the Villain model: the system undergoes two BKT type 
transitions with a massless U{1) symmetric phase at intermediate temperatures. At some 
particular "tricritical" value of h the massless phase shrinks to a point and it comes together 
with the first order transition line. This tricritical point is self-dual and is described by a 
conformal Z^ invariant parafermionic theory with the central charge c = 2{N — 1)/ (A^ + 2) 
introduced in In this type of models therefore generically one expects either the first 
order transition or a pair of BKT transitions with the massless phase in between. The 
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tricritical behaviour is special and requires fine tuning of tlie parameters. Tliis is indeed 
also the prevailing general expectation for the order of the transition in 2+1 dimensional 
gauge theories at large N: either first order or Villain type U{1) invariant behaviour. 

In fact we find our model in a completely different situation. The transition is not 
first order, and there is no U{1) invariant massless phase. We stress that within the RG 
flow eq.(|49|) the IR fixed point eq.(^) has two attractive directions. This means that it 
governs the IR behaviour of the points which lay on 2 - dimensional critical surface in 
the three dimensional parameter space, and is therefore generic. This by itself does not 



preclude that this fixed point is the same as the parafermionic Zn theory of |T^. If this 
is the case, it is quite interesting, since the point which appeared as "tricritical" from the 
point of view of usual spin models is in fact generic from the point of view of the 3D 
gauge theories. At present we can not prove that our critical point is described by the 
parafermionic theory but let us present some arguments supporting this conjecture. The 
point is that, as opposed to models considered in our Lagrangian eq. f^Of ) describes 
a theory of — 1 light fields. The theory of — 1 free massless fields have the UV 
central charge cuv = N — 1. However this CFT is deformed by the monopole and W - 
induced perturbations and flows to a different IR fixed point. However let us note that 
the central charge c = A^ — 1 is precisely the central charge of the SU {N)i WZNW model. 
The Ising (i.e. c = 1/2) model is the lowest among the minimal models with Virasoro (i.e. 
W2) symmetry. The highest model of this class is c = 1 model (one free field ) which is 
precisely SU{2)i WZNW model. When the c = 1 model was deformed by the monopole 
and W^-boson operators the central charge was reduced - and the resulting IR theory was 
Ising. 

Now, Z]sf parafermions with c = 2(A^ — 1)/A^ + 2 are the lowest minimal models with 
Wn symmetry - and the highest is SU{N)i (for more information about parafermions see 
for example |T3 and references therein) which can be described in terms of A^ — 1 massless 
fields. Thus if the theory in the UV describes N-1 massless fields and has symmetry, 
it is quite possible that result of the relevant (monopole+VT) deformation is a self-dual 
critical point. It is indeed known that the Zj\f parafermion theory is the self-dual model 
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with Wn symmetry. The fact that the central charge (and thus the effective number of 
degrees of freedom) is reduced in the process of the flow towards IR is of course in complete 
accord with Zamolodchikov's C-theorem. It is therefore possible that the IR fixed point 
that describes the universality class of the GG model is the conformal Z^v parafermion 
theory. 



Analysis of [T^, although admittedly incomplete also supports the expectation that 
we do not have Villain picture. In fact it is the presence of the large number of fields that 
drives our theory away from the Villain behaviour as we will now explain. 

4.3 Why not Villain? 

The RG equations eq.(^) were derived for the situation where all W bosons have equal 
masses (all fugacities fia are equal). One can wonder what happens if this is not so. In 
particular imagine an extreme situation, where some W bosons are light relative to the 
others, so that large monopole fugacities make all phase fields Xi (o^' components of ff) 
but one relatively heavy. In this case at zero temperature the theory seems to have only 
one light degree of freedom. This situation is as close as it can be to the spin systems 
with one phase field, and one may expect that in this region of parameter space the finite 
temperature behaviour will be similar to that in the Villain model. The appearance of the 
intermediate massless phase potentially has a natural place in our model. It could occur if 
the temperature at which the monopoles become irrelevant is lower than the temperature 
at which charges become relevant, Tbkt < T^^c- Then between these two temperatures 
the theory in the infrared is the theory of massless photons. Indeed, consider a simple 
invariant theory of one phase field 

C-eff = + eexp(2iV0) (51) 

We normalised the kinetic term so that for N = 2 the model reduces to the Polyakov 
effective theory for SU (2) GG theory. The BKT point in this theory is at 

'^BKT = ^ (52) 
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If the only vortices that are allowed have integer vorticity, the temperature at which they 
become relevant does not depend on and is 

T.C = (53) 

Thus for > 4 the "monopole binding" occurs prior to the "charge deconfinement" and 
there is an intermediate massless phase, bounded by two BKT transitions. 

Let us analyse in more detail how the model eq.(^DD behaves when one photon is much 
lighter than the others. The simplest case is SU{3) eg. (|3TD . Let us take Wi to be lighter 
than W2 and W3. This means that in eq. (|3lD we have ^1 ^ ^2, ^3- To minimise the first 
term in the potential, dynamically the difference of the phases of the two vortex fields must 
be constant. Thus on the low energy states we have 

Vi = V2 (54) 

With this identification we indeed get the theory of one phase field. However the coefficient 
of the kinetic term is "renormalised" due to the off diagonal form of Aij. In this case we 
find 

C,ff = ^i^,xr + ^exp{^3x) (55) 

This reduction procedure is easily extended to any A^. One can always choose appropriate 
W"s to be light, so that at low energy all vortex fields become equal 

Vi = Vj (56) 

The effective theory then is 

^eff = ^^iQ~2T^\ ^xr + ^eM^Nx) (57) 

Interestingly, the coefficient of the kinetic term of the only remaining field is of order A^^, 
which is the number of degrees of freedom in the underlying Yang-Mills theory. Thus the 
first thing to note is that the BKT temperature does not decrease as suggested by eq.(^), 
but rather increases with A^ as 



-BKT 



7fN 
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(5^ 



so that at — s> oo its value is twice that of = 2. 

To calculate T^c we should look at the terms that contain dual fields in eg . (|40|) . The 
structure of the phases in these terms is exactly the same as the structure of the phases 
in the monopole induced term. Thus clearly taking all Xi (and therefore Xi) to be equal 
some of these phases will vanish, while others will give the only surviving % field with the 
coefficient A^. Thus the charge terms reduce to 

/^exp(«^A^x)- (59) 



We then easily get 

g'N 



-NC 



(60) 



167r(A^- 1) 

So Tjvc decreases with A^. Perhaps surprisingly, we therefore find that as A^ becomes 
larger the two temperatures never cross, and in fact the difference between them grows. 
Nevertheless the temperature at which the scaling dimensions of the two operators are 

2 

equal always stays equal to the geometrical mean of the two temperatures in exact 
agreement with the analysis in the full theory eq. (|40|) . 

Why does this happen? If we were to allow only the vortices that preserve the con- 
dition Xi = Xji the only perturbations involving the dual fields would be of the form 

2 

fiexp{i^p^N{N — l)x)- This indeed would lead to much higher T/vc so that for A^ > 4 
the Tbkt and T/vc would cross. However the Lagrangian eq.(|40|) contains perturbations 
which create vorticity of a single phase field Xi, and thus effectively violate the equality of 
all phases. Another way of looking at it is to think of the field x in 6q.(0) as the average 
field X = Y,i!=i^Xi/^ ~ 1- The perturbations in eq.(|iO|) then induce fractional vorticity 
2tt/{N — 1). The corresponding operators are more relevant than those with vorticity one 
and thus the temperature T/vc is lower than one would naively expect. This effect is ob- 
viously due to the presence of the A^ — 1 independent fields all of which can be excited 
independently. Thus even though at low temperature the effective theory had only one 
light field, all fields are important in the transition region. 
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The preceding discussion is of course only illustrative, since it neglects the ef- 
fects of the lightest W bosons. Those light bosons lead to large monopole fugacity 
^ = exp{— 47rMi4//5f^}, which has an effect of freezing some of the phases of the vortex 
fields. However at finite temperature it is these same W bosons which are produced more 
copiously than the others due to their relatively large fugacity fi = exp{—Mw/T}. The 
appearance of these W bosons however tends to disorder precisely the same phase fields 
which are frozen by the corresponding monopole term by imposing non- vanishing vorticity 
on them. Thus the behaviour of the theory at criticality will be strongly affected by the 
presence of these particles and can not be directly deduced from the effective theory of 
only one scalar field, even allowing for fractional vorticity. 

It is interesting to note, that if we go high enough above the critical temperature 
where the monopole terms are irrelevant and can be neglected, the theory is described 
again quite well in terms of one light field. In this regime the large fugacity of light W^s 
leads to dynamical constraint Xi = Xj we have the theory of one light dual field. 

5 Relating to pure Yang-Mills. 

Although our analysis is not directly relevant to pure Yang Mills theory, it can be cast in 
the form which suggests that the relation exists and indeed may be closer than apparent 
at the first glance. 

The high energy phase of the Yang-Mills theory is indeed customarily described in 
terms of — 1 light fields. Those are the phases associated with the eigenvalues of the 
Polyakov loop, P [^. Since P is a special unitary matrix, it has — 1 independent 
eigenvalues. In fact these phases - the components of scalar potential Aq, are directly 
related to the dual fields rji of eq.(^OD |T5|. The dual fields fji appear in the last term of 



eq . (^Of ) . This term is nothing but the free energy of the charged particles W. This free 
energy is usually expressed in terms of P. In the regime where the Higgs expectation 
value is large and W^s are heavy, the only light components of the vector potential are 
the diagonal ones. Hence in this regime the Polyakov loop is naturally diagonal. The free 



19 



energy of a charged particle with the set of Abehan charges a is then given by the product 
of the appropriate eigenvalues of P. Comparing this with the last term of eq. (^0|) we have 

^0 = ^^^V (61) 

where Aq is the phase of the i—th eigenvalue of P. Remembering the the following relations 
between the roots and the weights of SU (N) 

N 1 

aij = w,-w,, J2 = 0^"' (62) 
i=i 2 

we can rewrite the effective Lagrangian in the hot phase (where the monopole terms are 
irrelevant) as 

rp N 

- HMf + Y.^^^^^os{A^ - A]) (63) 
9 i=i ij 

The phases exp i{AQ — A^} are eigenvalues of P, where V is the Polyakov loop in the adjoint 
representation. Eq.(pB|) is therefore of the form similar to the "effective action" discussed 
in the framework of hot QCD. Thus the "effective potential" in our case is given by a linear 
combination of the eigenvalues of the adjoint Polyakov loop. In fact, at the fixed point 
where all the fugacities are equal interestingly enough the potential term generated by W 
- can be written simply as 

/iTrP . (64) 

In the hot Yang Mills theory on the other hand the effective potential is given by the 
Bernoulli polynomial [Q. The origin of this difference is of course the large mass of 
W bosons in the GG model. The partition function of a heavy charged particle is well 
approximated by the Polyakov loop. Our derivation corresponds to the leading term in 
the low temperature expansion (expansion in powers of the Boltzmann factor) which in 
the GG model is valid even far above the critical temperature. In pure Yang Mills on the 
other hand the "charged particles" - gluons, are massless. As a result the particles are 
relativistic and their partition function is not given by the Polyakov loop. Also the low 
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temperature as such does not exist, and the standard perturbative calculation corresponds 
to the genuine high temperature expansion. 

Nevertheless it is interesting to observe, that some quantities calculated in the GG 
model behave in a way very similar to that in QCD. In particular consider the ratio of the 
longest correlation length in the sectors with total vorticity k, l/m^ to that of vorticity 
1. By total vorticity we mean the quantum number with respect to the magnetic 
symmetry. This correlation length can be extracted from the correlation functions of 
products of k vortex operators < Vi^...V^j. >. In general this calculation is quite laborious 
since the different vortex operators are not degenerate. However they do become degenerate 
on the trajectory leading to the fixed point, where all the fugacities are equal. As explained 
in |T5| at high temperature the inverse correlation length in the vortex channel is given 
by the "wall tension" of the Zn domain wall - solution of the equations of motion for the 
fields Aq with boundary conditions 

exp{iAl{x)} ^x^-oo 1, exp{iAl{x)} -^x^oo exY){i2-nk/N} (65) 

where x is the coordinate transverse to the "wall" . In the pure Yang Mills theory the result 



of this calculation is 21 



nik k{N-k) 



(66) 



mi N -I 

The equations of motion for the Lagrangian eq.(^) are (we take all variables to depend 
only on one coordinate) 
IT d"^ 

— ^(^0 - A^) + E ^^^^ sinK - ^^o] - E m sinK - ^^o] = (67) 

We are unable to solve these equations in the general case. However in two special cases 
they are easy to analyse. Consider first the case discussed in the previous subsection, 
when only one of the fields Aq is light. Then obviously on the solution we must have 
Aq = Aq = A, i,j,= 1,--,N — 1. Since A^ are phases of the eigenvalues of the special 
unitary matrix, the last component must then be A^ = (1 — N)A. These relations must 
hold for the solution with any k, including k = 1. Then A satisfies the equation 

2TN (f 

A + /xsin[A^A] = (68) 
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with fi = J2j^N f^Nj- This is the equation for one scalar field with potential cos[A^yl]. In 
this case clearly as long as k < N — k, the solution for 7^ 1 is just the set of k well 
separated solutions for /c = 1. When k > N — k, the same boundary condition eq.(^) can 
be satisfied by having N — k walls. Thus the tension of the /c-fold wall is 

rrik = min{k, N — k}mi (69) 



The other simple case is when all the fugacities are degenerate. Then following ||21 
we can try the following ansatz for solution 

Ai = kA, i = l,...,N-k (70) 
Ai={k-N)A, i = N -k + l,...,N (71) 

The resulting equation for A is 

2T (f 

A + fism[NA]=0 (72) 



This does not depend on k. The tension for such a solution scales as does the kinetic 



termpTI as ^(A^ — k). Thus the wall tension and the inverse correlation length in the 
channel with vorticity k scales like in hot Yang Mills theory according to eq. (|6^) . 

Thus even though generically the ratio mk/mi in the GG model is not universal, and 
depends on the details of the masses of the VT-bosons, close to criticality it follows exactly 
the same simple formula as in hot QCD. 

We can analyse in precisely the same way the behaviour of the ratios of the string 
tensions of fc-strings below the transition temperature. Due to the self duality of the fixed 
point, the effective Lagrangian in terms of the phases of the vortex operators Xi is identical 
to the Lagrangian for Aq with the substitution jj (, The tension of the 

confining string is then calculated as the tension of the domain wall separating vacua with 
different values of Xi 0- We thus find that the ratios of the string tensions also follow 
the relation eq.(^). In fact this scaling relation is commonly known under the name of 
"Casimir scaling" and is observed to hold for the ratios of the string tensions in pure 



Yang-Mills theory at low temperature p2[ in both four and three dimensions. 
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6 Conclusions 



An interesting feature of our result is that the critical temperature in the SU{N) theory 
at large is proportional to the coupling g"^ and not to 't Hooft coupling A = g'^N. Thus 
at large the critical temperature approaches zero. The physical reason for this is easy 
to understand. At large A^ and fixed A the Higgs VEV should also scale with A^ in such 
a way that the mass of W bosons remains fixed. The monopole action then grows as A^ 
and the photons get progressively lighter (exponentially with N)^. Thus the thickness of 
the confining string grows and the density of W bosons needed to restore the symmetry 
becomes smaller and smaller. 

More importantly, our main conclusion is that the deconfining transition in the SU{N) 
GG model is second order and the universality class is determined by the infrared fixed 
point eq. (|50D . This point is Zn symmetric and self dual. We have given some arguments 



supporting the possibility that the fixed point theory is the Zjy parafermionic model |18 



although we were not able to prove this explicitly. We can however definitely exclude 
Potts and Villain universality classes. In this context we also note that the ratios of the 
"wall tensions" calculated in the previous section (eg. (|66D) for A^ > 3 are different from the 
corresponding ratios in Villain model (which follow eg . (|69|) as well as in Potts model (where 
all the tensions are egual rrik = mi)^. This again tentatively supports our expectation that 
the universality class of the GG model is different. 

To answer this guestion one should study (numerically or analytically) the class of Zn 
invariant spin systems which has not been studied so far. The Lagrangian of the relevant 
model can be taken as eg.(^). This is an explicit Lagrangian of A^ — 1 interacting phase 
fields which can be easily discretized to define a lattice Zn invariant spin system. Hopefully 

the Wn symmetry of the SU{N)i WZNW model can be of help here too. 

^This is analogous to the situation in QCD where the instantons become less relevant at large N and the 
f]' meson becomes massless. The major difference is of course that while the rj' mass in QCD decreases as 
the photon masses in GG model decrease exponentially. This difference is due to the non diluteness 
of the instanton gas in QCD as opposed to diluteness of the monopole gas in the GG model. 

^Technically speaking the calculation of the previous section is valid only far enough from criticality, so 
that the monopole terms could be neglected. We believe however that due to the self-duality of the fixed 
point the same behaviour will also survive in the critical region. 
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Interestingly, contrary to naive universality arguments the transition is neither first 
order as in the A^-state Potts model, nor in the U{1) universality class as in the Villain 
model. We believe the reason is precisely the large number of light fields present in the 
theory. It is well known that oftentimes the symmetry alone does not fix the universality 
class of the transition, the number of light fields being the other important element. 

An interesting question is of course what happens in the pure Yang-Mills theory. The 
global symmetry associated with the phase transition is still Z]y [Q. The crucial question is 
what is the number of light degrees of freedom. We think there is some grounds to believe 
that the description presented in this paper is relevant in this case too. 

As discussed in the previous section there is direct correspondence in the hot phase 
between the light fields in the GG model and in the pure Yang-Mills theory. Again, the 
usual lore is that the behaviour of these same fields Aq at critical temperature determine the 
universality class of the transition. Moreover, the ratios of the vortex correlation lengths 
as well as string tensions close to criticality in the GG model seems to be similar to pure 
Yang Mills theory. This point of view would then fit with the proposition that the critical 
behaviour of the pure Yang-Mills theory is the same as that of the SU{N) GG model. Of 
course, universality arguments can never exclude the possibility of first order transition 
which can be forced upon the system by a heavy sector It would be interesting to 

investigate this question numerically by lattice gauge theory methods. 
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